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I. INTRODUCTION 


A most critical portion of flight, even in an ideal 
Operating environment, is the landing of an airplane at a 
desired landing site. Inclement weather at the landing 
Site, with the corresponding reduction in visibility, fur- 
ther complicates this portion of the flight operation to 
the point that landing may become an impossibility due to 
visibility restrictions. The increased demand, both com- 
mercially and militarily, for flight operations under all- 
weather conditions has prompted a great deal of research 
to develop an "Automatic All-Weather Landing System" capa- 
ble of automatically controlling an airplane during the 
landing portion of a flight. The problems encountered in 
such investigations, although not considered insurmount- 
able, are vast and complicated mainly because the flight 
dynamics of the airplane and the surrounding environment 
are themselves complicated. This investigation makes some 
simplifying assumptions to reduce this complexity and de- 
rives the necessary controls for such an all-weather land- 
ing system, using optimal control techniques. 

The problem and the simplifying assumptions are de- 
scribed in Chapter II. In Chapter III a general model for 
an airplane in the landing mode is developed, and the spe- 
cific model used in this investigation is presented. 
Chapter IV presents the optimal control theory employed 


and describes the method used to modify the original state 


variables to make the theory applicable. The problem 
specifications are introduced in Chapter V, the desired 
trajectory is formulated, and the significance of the per- 
formance measure in obtaining realistic results is dis- 
cussed. In Chapter VI the procedures used in the investi- 
gation are explained, and the cases investigated are pre- 
sented. The results are discussed in Chapter VII, and the 
conclusions and recommendations are presented in Chapter 


Vii. 


LO 


II. PROBLEM FORMULATION 


The basic descriptions and definitions of the landing 
portion of flight are presented in this chapter, along with 
the assumptions made to reduce the complexity of the pro- 


blem. 


A. STANDARD WEATHER CRITERIA 

The International Civil Aviation Organization (ICAO) 
has adopted the following set of minimum weather conditions 
for the automatic landing of an airplane, (Ref. 1). All 
distance specifications are given in meters, followed in 
parenthesis by approximate values in feet. 

Category I. Operation down to minima of 60 meters 
(200 ft.) decision height (altitude) and Runway Visibility 
Reading (RVR) of 800 meters (2600 ft.). 

Category II. Operation down to minima of 30 meters 
(100 ft.) decision height and a RVR of 400 meters (1200 ft.). 

Category III-A. Operation to and along the surface of 
the runway, with external visual reference during the final 
phase of the landing to a RVR minimum of 200 meters (700 
ites). 

Category III-B. Operation to and along the surface of 
the runway and taxiways, with visibility sufficient only 
for visual taxiing comparable to an RVR of about 50 meters 


WP5O Tt.) . 
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Category III-C. Operation to and along the surface of 
the runway and taxiway without external visual reference. 

The Federal Aviation Agency (FAA) has established the 
criteria for certification of any automatic all-weather 
landing system based on compliance with these categories of 
automatic landing weather minima. Presently no automatic 
systems that comply with any Category III minima have been 
certified. However, experimental tests have been success- 
fully conducted (Category III-B) with a system in a C-14l 
airplane (Ref. 2). Several automatic control systems have 
been certified for Category II operation, provided that the 
landing site has the prerequisite ground equipment. The 
vast majority of automatic landings are restricted to 


Category I conditions. 


Bs GROUND EQUIPMENT AT AN ALL-WEATHER LANDING SITE 

At present, only one special piece of ground equipment, 
an Instrument Landing Systems (ILS), is required for land- 
ing at an all-weather landing facility. The ILS provides a 
radio beam to establish a nominal glide path to the land- 
ing site and several locator beacons, placed along the land- 
ing track to provide distance checks. One locator beacon 
is positioned to provide an indication of the nominal deci- 
Sion altitude point for the facility. Both azimuth amd 
elevation information are provided by the ILS. Category I 
and Category II automatic landing operations are presently 


using this ground equipment with the latter, of necessity, 


eZ 


requiring more sophisticated ILS equipment to provide the 
increased accuracy needed in the radio beam. Figure l 
shows present requirements for the beam accuracy of ILS 
Category-II ground equipment. Several major U.S. air- 
fields have this equipment installed and operational. 
Because of beam distortions due to the electromagnetic dis- 
turbances, antenna installation environments, etc., the ILS 
equipment is not accurate below an altitude of about 60 
feet. Therefore, any automatic control of the airplane 
beyond Category II minima will require either additional 
ground equipment or an automatic system in the airplane it- 


self to effect the final phases of the landing. 





r — 24 ft 
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Figure 1 
Category II ILS Beam Accuracy Requirements 
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Cz PHASES OF ALL-WEATHER LANDING 

The automatic landing of an airplane under all-weather 
conditions can be divided into five phases: 

Le her App rOacCn = MeertCeuume mac’ 

The Approach-Intercept Phase is defined as that 
segment of the landing during which the configuration of 
the airplane is physically changed in anticipation of the 
landing and the airplane is directed toward the ILS beam 
which serves the landing site. 

2 Lem C apie — cack nace 

The Capture-Track Phase is defined as that segment of 
the landing during which the airplane is directed toward 
the landing point by an automatic flight control system 
which uses the ILS beam as its reference. 

3. The Flare Phase 

The Flare Phase is defined as that segment of the 
landing during which the airplane attitude is changed in 
preparation for the touchdown. 

4, The Land Phase 

The Land Phase is defined as the actual touchdown 

of the airplane on the landing surface. 
5. The Roll-Out Phase 

The Roll-Out Phase is defined as the climax of 

the landing during which the airplane is slowed to a stop 


or, equivalently, to a safe taxi speed. 
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Figure 2 depicts the phases of landing in both azimuth 
and pitch. Included, for clarity, is the decision altitude 
(height), which is defined as the altitude at which a deci- 
sion must be made as to whether it is feasible to continue 
the landing. As an example, the decision altitude for 
Category II landings is 100 feet above ground level. 

Automatic control of the airplane during the first 
three phases of the landing is presently standard operating 
procedure for both commercial and military flight opera- 
tions. However, lack of certified automatic flight control 
systems to accomplish the remaining phases of the landing 
dictates that the decision to continue the landing beyond 
the decision altitude be based on the premise that the 
pilot has visual contact with the landing site at the de- 
cision altitude, so that the airplane can be controlled 


manually, if necessary. 


D. ASSUMPTIONS MADE FOR THE LANDING PROBLEM IN THIS STUDY 
As stated in the opening remarks of this section, some 
simplifying assumptions are made in this study to reduce the 
complexity of the automatic landing problem. These assump- 
tions are listed below. 
Assumption 1. Only that portion of the landing from 
the decision altitude to the land phase will be considered. 


It is assumed that the airplane is automatically or manually 


lige: 


AZIMUTH 





our | Nd TRACK | INTERCEPT | 


FIGURE 2 


PHASES OF 
AUTOMATIC LANDING 


controlled, using the ILS as a reference, to the locator 
beacon which defines the nominal decision altitude point -- 
taken to be 100 feet above ground level. The roll-out 
phase is also assumed to be automatically or manually con- 
exolled. 

Assumption 2. The airplane is physically located in 
Space within the prescribed Category-II ILS "window" at the 
nominal decision altitude point and is in Equilibrium 
Flight (see definition in Chapter III) at that time. If 
these conditions are not met, it is assumed that the land- 
ing will be discontinued (airplane will be waved off). 

Assumption 3. Wind effects will not be considered. 
During landing, the airplane is subjected to both steady- 
state and gusty winds, the latter being of primary impor- 
tance, since the gusts are random in nature, Steady-state 
winds could be considered since their effect can be elimi- 
nated by a steady-state change in the airplane heading, 
whereas the random wind gusts could require a statistical 
handling of the problem. By neglecting both effects, the 
airplane velocity becomes equal to the ground velocity and 
the problem can be formulated in terms of time-to-go-to 
landing. 

Assumption 4. Only airplane motion in the vertical 
plane will be considered. Lateral motion of the airplane 
in the final phases of landing is primarily necessary to 


point the airplane in the direction of the runway just prior 
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to the actual touchdown. The airplane is normally "off" 


heading during the landing to counter steady-state cross- 


winds. Since the wind effects have been neglected, the 


lateral motion can be neglected. 
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Ifl. MODELING THE BLANT 


The first step in any control problem is the formula- 
tion of a realistic mathematical model to represent the 
dynamics of the plant to be controlled -- in this case, the 
airplane. In this chapter the formualtion of a general 
model applicable to any airplane is discussed, and the spe- 


cific model used in this investigation is presented. 


A. THE GENERAL MATHEMATICAL MODEL 

Aerodynamists have developed two sets of linear equa- 
tions which describe the dynamics of an airplane. These 
are referred to as the longitudinal or symmetric and the 
lateral or asymmetric equations of motion. Since, as pre- 
viously stated, this problem considers only the longitudi- 
nal motions of the airplane, further reference is limited 
to the longitudinal equations of motion only. Derivation 
of these equations is beyond the scope of this study, but 
complete and detailed derivations can be found in the li- 
terature. (Ref. 3, 4, and 5) 

Before presenting the equations of motion, it is sig- 
nificant, for clarification and reference purposes, to pre- 
sent several aerodynamic definitions and assumptions used 
in obtaining these equations. 

1. Aerodynamic Definitions 

a. Equilibrium sMieeheris *detmed "as unaccelerated 


flight. Flight is along a straight path during which the 


i 


linear velocity vector measured relative to fixed space 
is invariant and angular velocity is Zero. 

b. “Steady flight as defied as iliaghc) cucu 
which the linear velocity vector is invariant and angular 
velocity is constant. In this context, equilibrium, £Elicghe 
is also steady flight. 

c. Airplane coordinates, axes, and angles are 
defined in Figure 3, which shows the airplane in equili- 


Joneahibiy) az llavejales 


Equilibrium flight linear velocity along x axis 
a. = Equilibrium flight angle of attack 

oe = Equilibrium flight pitch angle 

9 Equilibrium flight glide angle 


Cig Instantaneous components are defined as the 


Summation of equilibrium flight components and correspond- 
ing perturbations caused by a disturbed flight condition. 
Figure 4 is a representation of the airplane in a disturbed 
flight condition and includes the perturbations of the two 


control components. 


V = Instantaneous linear velocity along x axis 

W = Instantaneous linear velocity along z axis 

Q = Instantaneous angular velocity along y axis 
A = Instantaneous angle of attack 

© = Instantaneous pitch angle 
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FIGURE 3 
AIRPLANE IN EQUILIBRIUM FLIGHT 
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FIGURE 4 
AIRPLANE IN DISTURBED FLIGHT 
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2. Assumptions Used in Obtaining Equations of Motion 


a. The airframe is a rigid body (no aeroelastic 
deflection of the airframe). 

b. The earth is fixed in space, and the earth's 
atmosphere is fixed with respect to the earth. 

c. The mass of the airplane remains constant 
during any particular dynamic analysis. 

d. The x-z plane (vertical plane) is a plane of 
symmetry. 

e. Disturbances from steady flight are small 
enough so that the products and squares of changes in ve- 
locities are negligible in comparison to the changes them- 
selves, and the disturbance angles are small enough so that 
the sines are equal to the angles in radians and the cosines 
are equal to one. Products of these angles are also approx- 
imately zero. Because these disturbances are small, the 
change in air density encountered by the airplane during 
any disturbance is considered to be zero. 

f. During steady flight conditions, the airplane 
is assumed to be flying with wings level and all components 
of velocity zero except for the linear velocity along the x 


axis. (Equilibrium flight) 
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g. The flow is quasi-steady. (The air flow pat- 
tern around the airframe instantaneously changes ina 
steady flow pattern as the airplane changes its orienta- 


tion with respect to its flight path.) 


3. The Longitudinal Equations of Motion 


Wa) = Kay it) + X a(t) - g@(t) + X w(t) 
i X W(€) i Xp 6,(t) + Kg S(t) 
w(t) = eget Vv ait) + Z a(t) + Z w(t) 
(3.1) 
+ Z(t) + 25 %'") + 25, on {*) 
Q(t) = Myv(t) + Myq(t) + M,w(t) + Mw(t) 


+ Me 6 (t) + M, 6. (t) 
6. 3 Sn at 
Where the subscripted X, Z, and M are the Dimensional Sta- 
bility Derivatives of the airplane which are parameters for 
a specific airplane in a specific flight regime (e.g., 
landing regime, subsonic regime, etc.). 


By substituting the relationships 


he 


w(t) vi &(t) 


(32) 


q(t) = 8(t) 
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the equations of motion can be rewritten in the more fa- 


miliar terms of v, @, and 9 which follow: 


v(t) = xX v(t) - g(t) + x ,8(t) + X q(t) 

+ XgB(t) + Xp GLE) + Xe Op (El 

Z ° Ae. Z 
&(t) = v(t) + (tb) + —* O(t) + = a(t) 

O Yo Yo 

Z Z 
Z 6 é (325 
; = &(t) + — 6 aaa 

Ve, We S we a 

Bit) = Mv(t) + M§(t) + Mya@(t) + Ma@(t) 


+Me 6 (t) + Ms 6 () 
= ah 

In a rigorous mathematical sense, these longitudinal 
equations of motion, which deal with perturbations from 
equilibrium flight, are applicable only to infinitesimal 
disturbances; however, aerodynamic experience has shown 
that quite accurate results can be obtained by applying 
these equations to disturbances of finite, non-zero 
magnitude. 

Since this problem deals with the landing approach, 
two additional parameters must be considered. One, the 


altitude of the airplane, is of primary importance, and an 
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equation relating actual airplane altitude to the flight 
dynamics is required. Such a relationship can be approx- 


imated by 
h(t) ¥ v T(t) (3.4) 
where h is defined as the instantaneous vertical distance 


above the ground of the aircraft wheels. By definition 


(see Figure 4) 


ey = y+ tt) 
(3) 4) 
y(t) = O(t) - a(t) 
and therefore, 
Ge ae Vv (8 (t) ~ (t)) (36) 
The second parameter is ground effect -- defined as the 


effect on the airplane dynamics as the airplane nears the 
vicinity of the stationary plane (the ground). This ef- 
fect on an airplane is non-linear and is dependent on the 
size and shape of the specific airplane. 

The formulation of a general mathematical model to 
describe the motions of any airplane in the landing ap- 
proach is now complete. All that remains is application 


to a specific airplane. 


ae THE SPECIFIC MATHEMATICAL MODEL 
The main problem associated with the use of the gener- 


al mathematical model is determination of the numerical 
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values for the Dimensional Stability Derivatives for the 
specific airplane. Because of the vast amount of experi- 
mental data available on the airplane, both from wind tun- 
nel studies and actual flight testing, the McDonald 
Douglas F-4J Phantom II jet fighter was selected. Appendix 
A contains representative data for this airplane in the 
landing configuration. 

As the F-4J Dimensional Stability Derivatives, aif Kaye 


lAgaer | ae ais 


r oO and Ms are zero and ground effects are neglected, 


a 
the longitudinal equations of motion and the height equa- 


tion can be simplified as follows: 


te X v(t) ee ae g(t) + Kg On (t) 
Ze i Ais oF 
a(t) a) 7 a(t) + 98(t) |= 6, (t) + opted 
(30) 
Q(t) =M v(t) + M,a@(t) + Mg a(t) ts M9 () * Mp 6,(¢) 
h(t) = vo (A (the- @(t)) tov 7, 
Defining 
V 
& 
x =| 86 (3.8) 
0 
h 
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and 


Ic 
lib 


op 


(3.9) 


the equations are manipulated into the state variable form 


x(t) S Arde) een Ec 


with the result 


ie “ioe “ee ° 0 
a5 aoy O ili O 
x(t) = xalate) 
0 0 0 i 0 
A471 342 or 0 
0 ao ap 3 0 0 
0 i 0 
Bo, = Po2 0 
+ suits) + 
0 0 0 
Ba, = Pg 0 
0 O Coy 


Za 


(Siako') 


uaalel: ) 


where 














ll Vv 52 O 
area aoe 
= ae Din = X% 
T 
a = ey b = 266 
Pal Vy. 21 V 
O O 
Z, 
_ 2a er: 
22 A 22 ve 
M32,, My § 
=A) a a, . V Day - Ms : V 
O e O 
MeZ Ma4 6 
= (Me ee b,. = = 
aA OL Vv A2 Vv 
O O 
243 a | ae C51 7 YoX% 


Appendix B contains the numerical values for these elements. 
Heretofore in the formulation, it was assumed that the 
two perturbation controls, elevator deflection and thrust, 
are instantaneously available to control the airplane. [In 
reality, both controls have associated dynamics, which are 
aerodynamically referred to as elevator actuator lag and 
lag between the thrust command and thrust. These effects 


can be approximated by use of the following differential 


equations: 
; 1 
b(t) =p — (8, (t) - 6 (t)) 
6 Cc 
: (3)2y 
: 1 
b(t) =o (6, (t) - 6,(t)) 
Op Cc 
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where the subscripted T's are the appropriate component 
time constants and 6. and bn are the elevator deflection 
and thrust commands sERPEGRITEIY: In this investigation 
it will be assumed that the controls are instantaneously 
available; i.e., the dynamics will be neglected. 

This completes the mathematical modeling of the F-4J 
airplane in the landing configuration. The model was test- 
ed on an IBM System/360 digital computer, using a Fourth- 
Order Runge-Kutta routine to solve the differential equa- 


tions. The responses resulting from various control in- 


puts were consistent with actual airplane responses. 
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IV. OPTIMAL CONTROL THEORY 


This chapter will not attempt to present the entire 
field of optimal control theory, but rather will assume 
that the reader is familiar with it; hence, only those 
techniques pertinent to this investigation will be dis- 
cussed. Some definitions will first be presented, follow- 
ed by mathematical techniques employed to obtain the op- 


timal control. 


A. OPTIMAL CONTROL DEFINITIONS 


lee ~COontro Tamns tony 


The history of control input values during the 
time interval L GE in is denoted by u. 
2.) obaltem mea Cemenmy 
The history of state values during the time in- 
terval hee is denoted by x. 
3. Admissible Control 
U denotes the set of all control histories which 
Satisfy the physical control constraints during the time 
interval [t ,t,J. 
4. Admissible Trajectory 
X denotes the set of all state trajectories which 
Satisfy the physical state variable constraints during the 


time interval Deemer! 
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5. The Performance Measure 
J denotes a scalar measure of the performance of 
a system when a control history is applied. 
6: Optirmaleconmtuet 
A control which causes the system to follow an 
admissible trajectory which minimizes the performance mea- 
sure is called an optimal control and is denoted by u”. 
7. Optimal Trajectory 
The admissible trajectory which results when an 
optimal control is applied to the system is called an op- 
timal trajectory and is denoted by x*. 
8. Tracking Problem 
A problem wherein the intent is to maintain the 
state trajectory as close as possible to a desired trajec- 
tory -- denoted by r -- in the interval [et is called 


a tracking problem. 


la8 THE LINEAR TRACKING PROBLEM 

Since the automatic landing problem as formulated in 
Chapter II is directed toward controlling the airplane in 
a desired manner during the final phases of a landing, and 
the mathematical model of the airplane as presented in 
Chapter III is linear, the problem can be considered as a 


linear tracking problem. 


onl 


It has been shown (Ref. 6 and 7), where complete deri- 


vations can be found, that, given a set of state equations 


K(t) = Ale) x ey) eee CUE) (ao 


where 


x(t) is the state vector of dimension n 
u(t) is the control vector of dimension m 
A(t) is a nxn matrix 

B(t) is a nxm matrix 


and the performance measure 


2 
a= 5 llx(t,) - x(t,) in 


e 
Fe 2 2 
+2] Ti xcey - ce? + dace? Jae 
t Q(t) Sly 


(4.2) 


where 
loc 2 L Jorn 
X(T) — ECM lo ¢gy = [XOT) ~ eT) | Q(T) X(T) - ECT) | 


the final time ty is fixed, x(t,) is free, u(t) 1S uncon= 
strained, H and Q(t) are real symmetric positive semi- 
definite matrices, and R(t) is a real symmetric positive 


definite matrix, that the optimal control exists and is 


unique. The optimal control is given by 


t) = F(t) x*(t) + g(t) (4.3) 
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where F(t) 1s the mxn matrix of feedback gains, and g(t) 
is the mxl command signal vector which is dependent on 
the system parameters and the desired state trajectory. 
Figure 5 is a block diagram of the plant and optimal con- 
teoller. 


F(t) and g(t) are given by 


F(t) - 


= R( Be BE ee act) (4.4) 


weil Be 


g(t)  - R(t) (t) s(t) (4.5) 


where K(t) is the solution of the Riccati-type matrix dif- 


ferential equation 


K(t) = - K(t)A(t) - A’(t)K(t) - Q(t) 
(4.6) 
+ K(t)B(t)R - (tyB. (t) K(t) 
with boundary conditions 
K(t,) = 3 (4.7) 


and s(t) is the solution of the linear vector-differential 


equation 
&(t) = - Cal(t)-K(t)B(t)R-"(t)B (t) Js (t) +0(t) x(t) 
(4.8) 
with boundary conditions 
s(t,;) = - Hx(t,) (4.9) 


i 


i 


FIGURE 5S 
BLOCK DIAGRAM OF PLANT AND 


CONTROLLER 
IN A LINEAR TRACKING PROBLEM 


Since the boundary conditions at the final time ty are 
known, these two matrix differential equations must be in- 
tegrated backwards in time from t to to: 


It was assumed that all states of the system can be 


measured, 


OF MODIFICATIONS TO SPECIFY STATE PLANT 
The previous discussion indicated that the optimal con- 
trol in a linear tracking problem exists and can be obtain- 


ed, assuming that the plant is of the form 
x(t) = A(t)x(t) + B(t)u(t) (4.10) 


Since the plant model, as presented in Chapter III, is in 


the form 
x(t) + Ax(t) + Bu(t) +c (4. ap) 


some modifications are necessary in order to adapt the 
specific problem to the theory presented. 

This modification can be accomplished by observing 
that the c vector contains only one non-zero element, 1.e., 
the differential equation describing the altitude of the 


airplane 


h(t) = Vo% + Vo (Ot) - a(t)) (4.12) 


contains the constant term va 


3 US) 


By defining 
h(t) # h(t) - h(t) (4.13) 
where 


fy UBS h(0) + vy (4.14) 


is the nominal equilibrium altitude of the airplane during 


the interval ie Cade It follows that 


h(t) = voy (4.15) 
pCR) STs) hf) (4.16) 
A(t) = Ut Vv, (8 (t) =O.) ) y= i iam (4.17) 
h(t) = v,(@(t) - a(t) (4.18) 


So h(t), as defined, is the perturbation of the altitude of 
the airplane about the nominal equilibrium altitude. 

With this modification, a new set of states, hereafter 
referred to as the revised states (to distinguish them from 


the actual states), can be defined as 


V O 
OL O 
x = 6 = x - O (4.19 
A 0 
h be 
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and the revised state equations become 
K(t) = AX(t) + Bu(t) (4.20) 


where the elements of A and B remain identical to those pre- 
viously formulated in Chapter III, equation (3.11). fMThis 
revised plant is now in the proper form for application of 


the optimal control theory. The solution then becomes 
a’ (te) = E(t)ege) + gic) (4.21) 


from which the behavior of the actual states can be deter- 
mined by substituting the definition of x(t). This substi- 


tution produces the optimal control 


u*(t) = F(t)x*(t) - F(t) | 0 + g(t) (4.22) 
; | 


h(t) 


By this simple mathematical manipulation, the problem can 
be solved by applying optimal control theory to the revised 
state equations and then expressing the results in terms of 


the actual states. 


oF 


V. SPECIFICATIONS, DESIRED TRAJECTORIES, 


PERFORMANCE MEASURE 


In this chapter the applicable specifications are in- 
troduced and a set of desired state trajectories for the 
problem are derived. Many of the results presented are 
based on the author's experience as a pilot operating the 
F-4 and other jet airplanes. The final section of this 
chapter discusses the significance of the performance meas- 


ure in obtaining realistic results. 


A. PROBLEM SPECIFICATIONS 

A successful automatic landing requires that certain 
conditions relative to the airplane and its environment be 
Satisfied. These conditions, which assume the role of spe- 
cifications for the problem, are formulated in terms of per- 
formance requirements and constraints on the system states 
and controls, together with specifications regarding related 
parameters. The following specifications are of primary 
importance to the problem under consideration. 

ee SC LOGEIEY. 

The instantaneous velocity of the airplane during 
the landing portion of flight must remain above l1.l Yani 
or V(t) 2 214.5 ft/sec. The upper velocity limit depends 
on the structural limit with landing gear and flaps extend- 


ed; flight experience indicates that a reasonable upper 
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Lanna we 1. Tt Vv, Ox V(t) s 245 ft/sec. In terms of the 
specified state, v(t), which is the perturbation velocity 


about equilibrium flight, these limitations are 


-8.5 ft/sec < v(t) < 22 ft/sec , Ee Labarat — 


E 


2. AnoglesvoeeAecack 


The angle of attack must remain below 0.9 A 
stall 
or A(t) < 0.44 radians. No stringent lower limit exists 
but, as in the velocity case, experience indicates that a 
reasonable minimum value is approximately 0.75%, or 


A(t) 20.25 radians. As a result, the perturbation angle 


of attack limits were established as 


lace) | < 0.11 rad , celta, tees « 


ii 


3. PitcheAngle 
During the landing portion of flight, the instan- 


taneous pitch angle is closely related to the angle of 
attack but does not have stringent limitations on its ex- 
cursions. However, at the actual touchdown point, limita- 
tions do exist primarily to prevent the airplane from 
either landing nose wheel first (low pitch angle) or tail 
first (high pitch angle). For the F-4J airplane these lim- 
itations, in terms of the perturbation pitch angle from 


equilibrium flight conditions, were established as 


-0.20 rad < 6(t) = 0.25 rae, Gl steest — 


£ 


o 


4, Pitch-Angle Rate 
Flight experience has indicated that this para- 
meter should be held to a minimum for pilot comfort. (See 
discussion under Related Specifications, below.) Realistic 


limits for the pitch angle rate were established as 


| (t)| < 0.08 rad/sec , Gqkae, Ged ‘ 


3. Altitude 
In the problem formulation, Chapter II, it was 
assumed that the altitude of the airplane at the initial 


time was within the prescribed Category II ILS "window" or 
SBi fitemeS h(t.) = a2" Ft 


In terms of perturbations about equilibrium flight condi- 


tions, this limitation can be expressed as 
hac eae 


As the airplane approaches the actual touchdown point, how- 
ever, altitude excursions about desired conditions should 
satisfy more stringent limits. As a result, the following 
limitations were established for altitude perturbations 
about a desired altitude trajectory, in the intervals in- 
dicated. (See Section B, this chapter, for desired trajec- 
EOL vey) 

In the interval from the decision altitude to the 


flare initiation altitude 


| h(t) = ha ()| S tere eet eee 


40 


anid in the inté@rval from the tlanre™inwtiatiwonsalltttuge to 


actual touchdown 
ln(t) - h(t)! = Sft, ae Ly 


where bane (ots) is the desired revised altitude state trajec- 
tory, and ty is the time at which the flare phase begins. 
6. Elevator Deflection 
Physical limits for elevator travel exist, and for 
the F-4J airplane in equilibrium flight in a landing ap- 


proach, were established as 
-0.26 rad < 6,(t) < 0.22 rad, Be hee te. 


7. Thimust 
Physical limits on thrust available exist for the 
F-4J airplane and are dependent on the selection of power; 
i.e., MILITARY engine operation or AFTER BURNER (A/B) en- 
gine operation. Flight experience has shown that realistic 
limits for thrust perturbations about equilibrium flight 
conditions in the landing portion of flight can be estab- 


lished as 
mSentee) | < 3000 lbs, ealit Sree | : 


This range is well within the thrust-available range of 


the F-4J airplane without A/B selection. 


4l 


Sve ~ Rell ditede Specifications 


Several related performance requirements and con- 
straints are applicable to the problem and are presented 
in the following listing: 

ayer Theat 't uel rouchdown™ pont mus tebe" wi chives 
feet of the desired touchdown point on the landing surface. 
Since the problem was formulated in the time domain, this 
restriction was established by requiring that the actual 
touchdown occur within + 0.65 seconds of the desired touch- 
down time. 

b. At the actual touchdown time the airplane must 
have no tendency to float, or, in other words, the airplane 
should have a positive rate of descent. In addition, the 
rate of descent at touchdown must be within the structural 
Sink rate limits for the airplane (Appendix A). Experience 
indicates that realistic rate of descent limits for the F-4J 


at the actual touchdown time are 
-9ft/sec < cee) < -3ft/sec 


c. The normal acceleration” GCEmEneC alrplanc.—— 
defined as a perturbation from the gravitational accelera- 


tion of equilibrium flight (1.0g) -- can be approximated by 


<j 


n(t) Sa Oe — @(t)) 


i} 
The airplane acceleration perpendicular to the hori- 
zontal reference plane. 
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and has physical limits based on the structural strength of 
the airplane. These limits are normally not attained in 

the landing portion of flight, due to the low velocity and 
control effectiveness in this regime. However, the func- 
tion n(t) gives an indication of the smoothness of the land- 
ing, and its effect can actually be felt by the pilot during 
the landing: therefore, to maintain the normal accelera- 
tion as closely as possible to the equilibrium flight condi- 


tion of 1.0g the following limitation was established 


| n(t)| < 0.29, Ee btamgieen 


f£ 
Bs DESIRED TRAJECTORY 

The problem has been formulated as a tracking problem; 
hence the desired state trajectory must be determined. This 
section presents a general desired trajectory appropriate 
for any airplane in the landing portion of flight, and in- 
dicates how these trajectories are applied to the F-4J 
landing problem. 

1. General Desired Trajectory 

Since the altitude of the airplane is of great 

importance during the landing, it is an obvious starting 
point in the formulation of a set of desired state trajec- 
tories. In Chapter II it was stated that this investiga- 
tion would consider only that portion of the landing from 
the decision altitude (100 feet) to the touchdown point. 


Referring to Figure 2, this portion of the landing includes 
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a section of the capture-track phase, the flare phase, and 
the land phase. In Chapter III it was shown that the pro- 
blem could be formulated in the time domain. These two 

considerations lead to a general altitude trajectory (spe- 


cified in the time domain) consisting of a constant descent 


from the decision altitude at t = ty = 0 to”-a Selected 
flare point at t = tie followed by a flare to touchdown at 
t = Cee 


In the optimal control theory presented in Chapter IV, 
one of the necessary requirements was that the final time 
(t_) be fixed. To alleviate the obvious disadvantage that 
this requirement imposes, it seemed appropriate to extend 
the final time beyond the actual touchdown time to ensure 
that aslanging.ocecurs prrenr. te t-- By so doing, an imag- 
inary flare plane below the actual landing plane was 
established. This implies that the airplane nominally will 
land before the final time specified. Figure 6 depicts 
this general concept in graphical form, where the capture- 


mace piase is the mum tcra-iam@eats the actual flare phase 


We 
is the interval [t,,t.), and the actual land phase is to. 
The imaginary flare phase is the interval [t,,t,) and the 
imaginary land phase is at the time ty. 
The rate of altitude change (rate of descent), h(t), 
must be considered in conjunction with the general altitude 


trajectory, even though it is not a state variable. Ideal- 


ly, during the capture-track phase LO, t,), the m@aterct 
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descent will be a constant, while during the flare phase, 
it will steadily decrease so that at the actual touchdown 
time (t.), the rate of descent will be within the estab- 


lished limits. Utilizing this consideration, the set of 


equations 
h(t) = h(t) + h(0) + ee er. 1) 
hate) = h(t) = v7, (5.2) 


which describe the desired altitude trajectory, h(t), and 
rate of descent, va a) for the interval [O,t,) were for- 
mulated. 

For the interval Lt tied , the desired altitude and 


rate of descent were defined as 


hg (t) 


L To (6 —t..) 
=?) l 1’. = _ 
h(O) + Vo%ot + ig (e 1 L,(t t,) 


(5.3) 


L (t-t ) 
ee Ce — -1 ) (5.4) 


ee) 
Since the revised altitude state was defined in Chapter IV 
to be 

h(t) = h(t) - h(t) (a5) 
equations (5.1) and (5.2) became 


Eee On EEO Es, (5.6) 
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Altitude 





O 
Figure 6 
General Desired Altitude Trajectory 
and 
_ L L, (t-t, ) 
ae ik Ulisse hee 

hy (t) S ~1) L, (t-t,) C5 a7) 
= L, (t-t,) 1 
hg (t) = Ly Ce = ; telt»t, : (S83) 


Equations (5.6) and (5.7) represent the desired revised 
altitude trajectomysinsihe si hiten a | [o,t,] in terms of per- 


turbations about the equilibrium altitude. 
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These general expressions for the desired revised rate 
of descent, hee) provide insight into the form of the 


desired angle of attack trajectory, Ose and pitch angle 


trajectory, Oa since from equation (4.18) 
h(t) =v (6(t) - a@(t)) . (5.9) 
Therefore, 


6 5 (t) SCA) 0, telO,t) (5. 16) 


and 


6 5 (t) = eet) y(t de) 


Ti, Leese) 
<aie 1 1 ils) , telt.,t 
Wes Ae: 


By definition (Chapter II), ® and @ are zero when the air- 
plane is in equilibrium flight, hence equation (5.10) is 
consistent. The problem remains to find equations for @ 
and 8 in the interval [t,,t,! which are consistent with 
the requirements of equation (5.11). Before proceeding, 
recall that the desired pitch-angle-rate trajectory, 9 
must also be considered. An acceptable expression for the 
desired pitch-rate trajectory is a linear function during 


the flare phase interval Ree given by 


6 4 (t) Senta CERI (anes 


1? 
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Integrating equation (5.12) and applying the desired bound- 


ary condition aty te= tiwei.e., 8 (t,) = 0 
L 
aes) 4 
OE) =) eee (5 1p) 


Substituting equation (5.13) intees.]Peeields 


L L L, (t-t 
- rn 2 _ 2 1 
a(t) = 5 (t t,) 7. (e 


) 
* =1) (5) 


as the expression for the desired angle of attack in the 
interval Be o Brcie 

The final state component trajectory which must be 
considered is the velocity trajectory. Experience from 
both an operational and a safety standpoint has establish- 
ed that the airspeed during a landing should remain es- 
sentially constant. This implies that Vge the desired 
velocity trajectory, should be zero during the entire in- 


terval of interest, hence 


vatt) = 0, ae ee 2 (5.15) 


2. Specific Desired Trajectories 


In the preceding section a general set of desired 
State trajectories -- applicable to any airplane -- were 
formulated. What remains is application to the F-4J air- 
plane so that the unknown constant terms can be evaluated. 
Since, in the interval [0,t,), the desired state trajec- 


tory is zero, the following discussion will be concerned 
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with the interval Leiatee where the equations of interest 
eC wl 5.3) (5.4), (5.7 ee (5 Sede (Saleen (50 teem a (5 59h4) 
-- hereafter referred to as the reference equations. 

The first consideration, however, was to determine a 
realistic time frame for the problem. The equilibrium air- 
speed, Vor of the F-4J, as given in Appendix A, is 223 
ft/sec. Assuming that the airplane is at the nominal de- 
cision altitude (100 feet) at t = 0 and that the airplane 
remained in equilibrium flight on a negative three~-degree 
glide slope, the actual time of touchdown from equation 
(4.14) is approximately 8.57 seconds. This assumes that 
there is no flare phase. Although the standard operating 
procedure in the U.S. Navy is to operate the F-4J airplane 
with no flare, this investigation will include a flare 
which commences at ty = 6 seconds; at this time the air- 
plane is at a nominal altitude of approximately 30 feet. 
It follows then that a reasonable time frame for the prob- 
lem is 10 seconds, and that the interval (ies piers is 4 sec- 
onds in duration. 

Next, values for the unknown constants in the refer- 
ence equations must be obtaitneas™*a@tmaary Sonametons aor 
the reference equations at t = ty = 6 seconds are known, 
and final-time boundary conditions can be selected to pro- 
vide realistic trajectories. However, recalling that the 


actual touchdown of the airplane occurs prior to the final 


time, and that the specifications refer to this actual 
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touchdown time, another set of boundary conditions at 

t= ty must be satisfied. The conditions at t = t. are not 
exactly given, but realistic values can be selected to meet 
the necessary See EL OmenOne As a result, values for the 
unknown constants Li, Lo, and L. of the reference equations 
were found by a trial-and-error method. Values of these 
constants were arbitrarily chosen to satisfy boundary con- 
ditions at t = te and then substituted into the reference 
equations to obtain the actual touchdown time, to, and the 
corresponding values of the parameters at t = to. This 


process was continued until values of Lie Lo, and L. were 


found, which resulted in realistic trajectories. The values 


Ly = 0.25 
L. = 4,68 
L. = 0.016 


which were obtained result in a desired actual touchdown 
time, (t.), of 9.3 seconds. The reference equations which 


apply to the interval [6,10] are then 


7 ARES . 25'(E=6) a _ . 

Dee) ee eee 1) 4.68 (t-6) 
( 55916) 

; Z AEH ESS) 

heii) n=) vilagimaty 462 (e = (.Sicullatyy) 
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Aa(t) = eA) Ga ai Be 68 (t-6) (5.18) 


0.25 

h(t) SC Came (5.19) 

@4(t) = 0.016 (t-6) (5.20) 

64(t) = 0.008(t-6) * (5.21) 

a(t) = 0.008(t-6)* - 4:85 (e°“°Mi.1) (5.22) 
0 


while in the interval [0,6) the applicable equations are 


h(t) = h(O) + voy = (5.23) 

h(t) = wey, (5.24) 

ey eee) = O(c) — ate) Semen 0 (5.25) 
and 

aml EE (i 3.0) a (Be 26) 


Table A summarizes the solutions of equations (5.16) through 
(5.22) at the flare initiation time (t,), the actual touch- 


down time (t.), and the final time (t.). 


Sel 


Table A 
Solutions to Reference Equations at Flare Initiation 
Time, Actual Touchdown Time, and Final Time 


Solutions at 


Parameter Equation No. t=6 sec t=9.3 sec t=10 sec 


h, (£t) 5.16 29.947 0 ~3.286 
h, (£t/sec) 5.17 -~11.676 -~5.667 -3.621 
h,(£t) 5 Ls 0 8.582 13.469 
h, (£t/sec) 5,19 0 6.009 8.055 
6 ,(rad/sec) 5.20 0 0.053 0.064 
6 , (rad) Beil 0 0.087 0.128 
a , (rad) Se 2 0 0.060 0.092 


Figure 7 presents a plot of the actual desired altitude 
trajectory with the equilibrium altitude trajectory includ- 
ed for comparison. Figure 8 illustrates the revised desired 
altitude trajectory Aas while Figure 9 depicts the three 


angular components of the desired state trajectory, ass 


04 and 5 and the velocity trajectory v.. 
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© 


time-sec 


Figure 7 


Desired Altitude Trajectory 


time-sec 


Figure 8 


Desired Revised Altitude Trajectory 
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time-sec 


Figure 9 


Desired Velocity, Angle of Attack, Pitch Angle, 
Pitch Rate Trajectory 
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and 


Ga PERFORMANCE MEASURE 
The mathematical form of the performance measure in- 


troduced in Chapter IV was 


2 
J =1/2 | x(t,~) - x(t,) | 
H 


ic 
fe 2 2 
+5 | [il xcey - xe) + lacey] at 
t Q(t) R(t) 


O 


(5.27) 


One of the advantages of this quadratic performance measure 
is that the elements of the H, Q(t) and R(t) matrices,which 
shall be called the weighting matrices, can be related to 
the design parameters of the system and chosen to satisfy 
the design objectives. If these weighting matrices meet 
the requirements established in Chapter IV, specifically 

H and Q(t) are positive semi-definite and R(t) is positive 
definite, the optimal control law can be found, and is 
unique for the selected set of weighting matrices. 

Previous investigations dealing with similar airplane 
landing problems (Refs. 8 and 9) have indicated that the 
weighting matrices must be diagonal matrices with positive 
diagonal entries. In other words, the weighting matrices 
must have non-zero weighting values specified for each 
state and control. As a result, in the initial trials it 


was assumed that the weighting matrices are diagonal and 


a 


time invariant over the interval of interest. Since the 
form of the weighting matrices was so selected, each dia- 
gonal element of the respective matrix could be related to 
a specific state variable or control, and therefore, theo- 
retically, could be assigned a proper value to meet estab- 
lished specifications. Unfortunately, this selection of 
weighting factors was not simply accomplished. In fact, 
it became obvious during the investigation that the state 
trajectory was very sensitive to changes in the weighting 
matrices. (See related discussion in Chapter VII.) 

: 

The selection of values for the weighting matrices was 
carried out by a trial-and-error method wherein a set of 
weighting matrices was selected, the optimal control law 
computed, the optimal trajectory calculated, and the re- 
sults compared with established specifications. In effect, 
the final weighting matrix selection was based on obtain- 
ing a realistic optimal trajectory which conforms to the 
specifications. 

Appendix B contains the numerical values used in the 


weighting matrices for this investigation. 
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VI. INVESTIGATION PROCEDURE 


Since the previous discussion has presented the pro- 
cedures used in modeling the plant, developing the desired 
trajectory, and related topics, only the procedure used to 
obtain the optimal control will be presented. The actual 
cases investigated will also be discussed in the last sec- 


tion of this chapter. 


A. PROCEDURE 

The following procedure was used to obtain the op- 
timal control: 

1. A set of representative values were selected for 
the diagonal elements of the weighting matrices (H, Q, and 
R) of the performance measure, equation (5.27). 

2. Equations (4.6) and (4.8) were simultaneously in- 
tegrated from t = ty to t = 0 to obtain values for K(t) 
and s(t) respectively. 

3. Using the results of step 2, equations (4.4) and 
(4.5) were solved to obtain values for F(t) and g(t) re- 
spectively. 

4, The state equations (3.11) were integrated from 


t=O tot=t using the results of step 3 to obtain 


‘has 
wy (t) of equation (4.22). 
5. The optimal control and optimal trajectory ob- 


tained from step 4 were observed to ascertain compliance 


5 


with problem specifications. If the results were not ad- 
missible or unrealistic, the values of the elements of the 
weighting matrices were changed and steps 2 through 4 were 
repeated. 

To implement the foregoing procedure, a computer 
program was written in Fortran IV for use in the IBM Sys- 
tem/360 digital computer. System/360 Scientific Subrou- 
tines were used in the program to accomplish the matrix 
algebra. An existing Fourth-Order Runge-Kutta subroutine 
developed at the Naval Postgraduate School was used for 
the required numerical integration. The computer program, 
including the integration subroutine, is presented after 


the appendices. 


Be CASES INVESTIGATED 

Two cases were investigated in the initial study. 
Case I assumes that the velocity of the airplane is con- 
stant and elevator deflection is the only control avail- 
able; Case II considers the plant and controls as pre- 
viously formulated. As a result, two different state mo- 
dels were used. 

1. Case I 

Since it is assumed that the velocity is constant 

and that elevator deflection is the only control available, 


the state model, equation (3.11) reduces to 


58 


Bae) = x(t) 


(6.1) 


Ze 


4l 


oylk 


where 


x = (6.2) 


and 


ue «= 6 Von) 
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By the same argument as was presented in Chapter IV, this 


state equation can be manipulated into the revised state 


eguation 
*(t) = AX(t) + Bu(t) (6.4) 
where 
CO 
- 6 
x = : (6375) 
6 
h 


and the elements of A and B remain identical to those of 
equation (6.1). 

2. Case II 

The state model, equation (3.11), and the revised 
state model, equation (4.20), as previously presented, were 
used in Case II. 

In order to evaluate the optimal trajectory, three dif- 
ferent sets of initial conditions were selected. The first 
assumes that the airplane is in an ideal flight condition, 
or, in terms of the revised state model, all initial state 
values are zero. The other two sets were selected to con- 
form to the worst, initial, flight condi t10nssthnat couldsoe 
reasonably expected. The following listing discusses 


these initial flight conditions. 
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a. HIGH AND FAST. This set of initial conditions 
assumes that the airplane was at the upper altitude limit 
of the ILS Category II "window" and that the airplane ve- 
locity is above Vor the equilibrium airspeed. In Case I, 
the second assumption is not, strictly speaking, applicable 
Since it is assumed that velocity was constant. However, 
the condition was simulated by assuming that the initial 
angle of attack, a(t), and pitch angle, B(t), were below 
their respective equilibrium values. This condition would 
result in actual flight if the airplane velocity were 
above the equilibrium flight velocity. It was further 
assumed that the initial pitch angle rate, Ot Dy was zero. 

b. LOW AND SLOW. This set of initial conditions 
assumes that the airplane is at the lower altitude limit of 
the ILS Category II "window" and that the airplane velocity 
is below the equilibrium airspeed. By the same reasoning as 
before, the second assumption was simulated for Case I by 
assuming that the initial angle of attack and pitch angle 
were above their respective equilibrium values. Again, 
the initial pitch angle rate was assumed to be zero. 

Table B presents a summary of the cases investigated 


and includes the initial values used for each case, 
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Table B 


Summary of Cases Investigated 


Bl Eval rli ght Applicable Initial State 
Case Gondisk1on State Equations Conditions 
IA IDEAL (6.1/6.4) (i; OmO aw) 45 
IB HIGH AND FAST ( SualGwA)ee(e0 Soee 2ovcemam 12m) 
Ic LOW AND SLOW (6.1,6.4) (0.03,0.03,0,88.0) 
IIA IDEAL (3.11,4.20) (0,0,0,0,0)7 
IIB HIGH AND FAST (Cevyeow20)MKSMo 08026 40k CoKoendan: 
IIC LOW AND SLOW (Gel ledw2o) (abeOn10K0Sn0N03 AOMReRO iE 
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VII. RESULTS 


A general result is first presented, followed by spec- 
ific results, in graphical form, for each of the cases in- 
vestigated. The final section of this chapter discusses a 


problem encountered during the investigation. 


A. GENERAL RESULTS 

The optimal control for an all-weather landing in the 
F-4J airplane was derived (see Specific Results below) 
using the techniques presented. Since some simplifying as- 
sumptions were made to reduce the complexity of the prob- 
lem, the results are not conclusive, but serve to demon- 
strate that the design of an automatic controller for the 
landing of an airplane is feasible by formulating the prob- 
lem as a tracking problem and applying optimal control 
theory. By systematically eliminating some of the assump- 
tions made in this study, a closer approximation to the 
actual airplane landing problem can be accomplished. 

Actual implementation of an automatic controller to 
provide the optimal control derived by these techniques 
would require that a digital computer be placed aboard the 
airplane which may impose an undersirable penalty. If so, 
Ene eptimal control still prewmdes "ene *coummsor the de- 
Sign of any sub-optimal controller and hence is of great 


value to the control engineer. 
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By SPECIFIC RESULTS 
Specific results for each of the cases investigated 
are presented in graphical form and are summarized below. 
1. The solution to equation (4.6), the K(t) values, 
are presented for: 
a. Case I in Figure 10 
b. Case II in Figure 21 
2. The solution to equation (4.8), the s(t) values, 
are presented for: 
a. Case I in Figure 11 
b. Case II in Figure 22 
3. The optimal control is, depicted for: 
a. Case IA in Figure 12 
b. Case IB in Figure 15 
¢. CaseygIC in Figure 18 
ad. Case IIA in Figure 23 
e. Case IIB in Figure 26 
f. Case IIC in Figure 29 
4. The optimal and desired altitude trajectories are 
presented for: 
a. Case IA in Figure 13 
b. Case IB in Figure 16 
Coe Casemic in agure: 19 
d. Case IIA in Figure 24 
e. Case IIB in Figure 27 


f. Case TIC in Frourers? 
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5. The components of the optimal trajectory, v, 
a, 9, and 6, together with the related airplane parameters, 
rate of descent and normal acceleration, are presented for: 

a. Case IA in Figure 14 

b. Case IB in Figure 17 

c. Case IC in Figure 20 

d. Case IIA in Figure 25 

e. Case IIB in Figure 28 

£. Case TIC in Figmre 31 

In all cases, the results were within established 

specification limits and were considered realistic for an 
airplane landing. The anticipatory nature of the optimal 
control was evident in the results, especially in Case IIA, 
where the optimal control produced a finite thrust pertur- 
bation at the initial time, although the airplane was as- 
sumed to be in an ideal flight condition. This apparent 
lack of continuity, along with the discontinuity encounter- 
6a) in the optimal control of Cases IB, IC, IPMS, ana ire at 
t = 0, was predictable, since it was assumed that the con- 
trols were instantaneously available. Development of a 
state model to include control dynamics, presented in 
Chapter III, will presumably eliminate this apparent dis- 


crepancy. 
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o PROBLEM AREA 

The importance of proper selection of the constant 
values for the elements of the diagonal weighting matrices 
of the performance measure is evident; namely, no unique 
values exist and hence proper selection is based on pro- 
ducing an admissible control and trajectory when the op- 
timization technique is applied. In this investigation, 
where many specifications had to be simultaneously satis- 
fied, this selection became paramount, especially when it 
became apparent that several components of the state tra- 
jectory were sensitive to changes in one element of the 
weighting matrices. As a result, the trial-and-error pro- 
cess required to obtain values for these weighting matri- 
ces was tedious and time-consuming. Further study is 
warranted toward the development of a scheme or method to 
aid the control engineer in the selection of these weight- 


ing matrices. 
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VIII. CONCLUSIONS 


For the results presented in Chapter VII, the follow- 
ing conclusions are stated: 

1. The feasibility of obtaining realistic results by 
formulating the all-weather landing problem as a linear 
tracking problem and applying optimal control techniques 
was demonstrated. 

2. Satisfactory landings, using a simple mathemati- 
cal model of the F-4J airplane, were accomplished in all 
cases investigated when the derived optimal control was 
applied. 

3. Obtaining values for the elements of the weighting 
matrices, H, Q, and R, to produce admissible and realistic 
results was difficult because the state trajectory was 
sensitive to changes in these matrices. The process was 
problematical, since a large number of design specifications 
had to be satisfied. 

Recommendations for logical extensions of this investi- 
gation follow: 

1. Formulate the problem to include the control dy- 
namics. 

2. Formulate the problem to include measurement noise 
and wind effects. 

3. Investigate application of sub-optimal control 


techniques to the problem. 
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4. Investigate the possibility of using the sensi- 
tivity of the state trajectory to changes in the weighting 
matrices as a basis for developing a method of obtaining 


representative values for the elements of these matrices, 
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APPENDIX A REPRESENTATIVE DATA FOR THE F-4J AIRPLANE 


IN THE LANDING CONFIGURATION 


Gross Weight a2, 500° Ib 
Equilibrium Velocity (vo) 223 ft/sec 
Stall Velocity (Vo¢311) 195 ft/sec 
Equilibrium Angle of Attack (a) 03337 rad 
Stall Angle of Attack (A stall 0.49 rad 
Limit Sink Rate at Touchdown 23 ft/sec 
Approximate Thrust Available 
MIL engine operation 5,000 lb 
A/B engine operation 12,000 lb 
Aerodynamic Ground Effect Negligible 
Control Time Constants | 
Elevator (T¢» ) 0.10 sec 
e 
Thrust (Tee } 0.15 sec 
e 
Dimensional Stability Derivatives 
X =-0.593x10 ‘sec? Za = 0.0 
Vv 
X = 0.0 Z6 =-0,132x107ft sec? 
q & 
X,= 0.107x10°ft sec™ Zp,= 0-289x10°°ft 1b™ 
Xa= 0.0 Wh = 0.450x10°-£t ‘sec 
Xp = 0.0 M_ =-0.430 sec + 
e qd 
X4,=0.950x10 "ft lb7+sec~? My =-0.130x10 sec 
Z., =-0.259 sec > Ms =-0.260 sec + 
7 -2 
5 Blak M§,,=-9-279 sec 
Z =-0.827x107£t sec? Ms = 0.0 
a T 
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APPENDIX B SUMMARY OF NUMERICAL VALUES 


STATE MODEL CONSTANT 
i 


= = -0.59300x10 > aq> = 
ee 0.10700x10° ins = 
ay, = -0.32172x10 Dan = 
a5, = -0.11628x10~* b5; = 
Aao5 = -0.37085 bo» = 
a,, = 0.75232x10~~ Dane 
ano = -0.12036x10 bao = 
ang = -0.69000 —lo ia 


DIAGONAL WEIGHTING MATRIX CONSTANTS 


CASE I 
aie BOS Sa 
Dog > Lol Ely = 
h,, = 1.0 G33 = 
By) caesiolee dis = 
R = 5.0 

CASE IT 

hj, = 5.0x10 > qi = 
iy = ola doo = 
Bee) SOOO das = 
me ee dag = 
Aes = 5.0x10 7 des = 
ri = 5.0 Loo = 


a2 


-0.22300x10°7 


0.22300x10° 
0.94667x10 ~~ 


-0.59193x10 72 


0.12979x107> 


-0.27746x10 


-0.33745x10 © 


OMG 6 1G 


1.0x107+ 
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